Maximal work extraction from quantum systems 
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Thermodynamics teaches that if a system initially off-equilibrium is coupled to work sources, 
the maximum work that it may yield is governed by its energy and entropy. For finite systems this 
bound is usually not reachable. The maximum extractable work compatible with quantum mechanics 
( "ergotropy" ) is derived and expressed in terms of the density matrix and the Hamiltonian. It is 
related to the property of majorization: more major states can provide more work. Scenarios of work 
extraction that contrast the thermodynamic intuition are discussed, e.g. a state with larger entropy 
than another may produce more work, while correlations may increase or reduce the ergotropy. 



The generality of the laws of thermodynamics for 
macroscopic bodies led to a long-lasting effort to derive 
them from microphysics 1,3,4 . This program is by now 
completed. In contrast, for finite systems, the applica- 
tion and even the formulation of the laws of thermody- 
namics are still the subject of studies^. The origin of 
qualitative differences between large and finite systems 
was recognized long a.go£&. A quantum system submitted 
to time-dependent external potentials that describe work 
sources undergoes a unitary transformation. During such 
an evolution, the density matrix has constant eigenval- 
ues, thus it cannot become gibbsian when starting from 
an arbitrary initial state. In contrast, most macroscopic 
systems have a thermodynamic behavior: they evolve 
close to a Gibbs state under the effect of slowly varying 
external potentials, keeping nearly constant entropjA2. 
This holds in general, although their evolution is Hamil- 
tonian at the microscopic level, owing to the large and 
smooth density of states of the spectrum^. Accordingly, 
the responses of finite and of infinite systems to exter- 
nal perturbations differ qualitatively. As an example, 
processes were studied which would be reversible in ther- 
modynamics, but entail a specific irreversibility due to 
the finite size of the involved systemifliii. 

The problem we treat here is well known in thermo- 
dynamics. It initiated its birth in 1824: What is the 
maximal amount of work that can be extracted from a 
system S by means of an external source of work acting 
cyclically in a thermally isolated process? And what are 
the criteria for comparing different states with respect to 
their work-providing ability? We find answers to both 
questions, and show that for finite systems they qualita- 
tively disagree with the standard thermodynamical ones. 

A specific case in which a finite system S behaves ther- 
modynamically in regard to its work production is well 
known. If S is initially prepared in a gibbsian state, it 
cannot produce work when coupled to a source of work 
through an external variable which returns to its initial 
valu o 4 ! 9 ! 12 . Such a state is thus called passive. We shall 
therefore be interested in off-equilibrium initial states. A 
scenario for preparing such states is to take S consisting of 
two non-interacting (or weakly interacting) subsystems, 
to set separately the two parts in thermal contact with 



heat baths at different temperatures, and then to decou- 
ple them from the baths. Subsequent coupling with a 
source of work may yield an amount of work that we 
wish to evaluate. 

The maximal work-extraction problem is thus posed in 
the following waji£. Consider a system S which can 
exchange work with external macroscopic sources. The 
evolution of its density operator p(t) is then generated 
by a Hamiltonian H(t) = H + V(t), where the time- 
dependence of V(t) accounts for work transfer. Following 
Refsi^, we call "cyclic" a process in which S, originally 
isolated, is coupled at the time t = to external sources 
of work, and decouples from them at the time r. Thus, 
the driving variables of the sources are cyclic, and the 
potential V(t) vanishes before t — and after r: V(0) = 
V(t) = 0. However, S need not return to its initial state 
at the time r. The system S is thermally isolated but may 
involve energy exchanges between its parts. The initial 
state p(0) = po and the Hamiltonian H being given, we 
look for the maximum work W that may be extracted 
from S for arbitrary V{t). According to the dynamics 

ih P =[H(t),p(t)}, (l) 

the work dW = tv[p(t) V{t) ] dt done on S isi 2 - 3 . 

W = tr[ P ( T )H] - E(p ), E(p Q ) =ti[p H] = E x . (2) 

Among all final states p(r) reached from po under the 
action of any potential V(t), we are therefore looking for 
the one with lowest final energy E[ = tr[p(r) H). 

The standard answeri^ relies on the idea that the final 
state p(r) is Gibbsian and that its von Neumann entropy 
S = — trplnp cannot decrease between the times and 
r. The maximum value of the work (J2J) is reached when 
the final state has the equilibrium form 

e -0H 

p(r)=p eq =—, Z = tre-? H , (3) 

with P > determined by the equality of the initial and 
final entropies: In Z — f3d\n Z/d(3 — — trpohipo- The 
largest amount of work — W extractable from S is thus 

W th = E(p )- TS(po) +T\nZ, (4) 
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the familiar difference of free energy between initial and 
final state, both evaluated with the final temperature T. 

Finite systems. The above derivation involves two ar- 
guments which call for some discussion. Following ther- 
modynamical intuition, we have first stated that the en- 
tropy cannot decrease. In fact, the von Neumann entropy 
S(p) remains constant during the evolution If S is 
macroscopic, it is the coarse-grained entropy which can 
increase. Anyhow, this point is harmless since the bound 
Q corresponds to constant-entropy processes of the to- 
tal system. We also implicitly assumed that S may be 
brought into an equilibrium state p cq by means of some 
evolution Q). For macroscopic systems this is usually 
allowed since dissipative processes within S may occur 
while it evolves under the influence of the coupling V(t). 

For finite systems, the sole action of V{t) is in general 
not sufficient to allow reaching at the time r a gibbsian 
state of the form J3J). Indeed, not only is the entropy S(p) 
conserved during the evolution generated by Q), but all 
the eigenvalues of p*&. In contrast to thermodynamic 
systems, finite systems keep memory of their initial state 
and do not involve any relaxation mechanism. One may 
therefore expect that the maximal amount of work W 
extracted from S is generally smaller than WW 

More precisely, the evolution Q of p is unitary so 
that p(r) = U po UK We look for the minimum of the 
final energy E[ ~ tiU poU^H over all unitary opera- 
tors U. We can parameterize the variations 5U of U as 
SU = XU , where X is an arbitrary infinitesimal antiher- 
mitian operator. Hence, we find SE{ — tr(XUpolI' H — 
Up WXH) = irX [p(r),H]. The stationarity of E { thus 
implies that p(r) should commute with H and have the 
same eigenvalues as po, a condition which replaces J2J). 
In the spectral resolutions of po and H, 

P° = ^2 r j\ r i)( r j\> H = ^2e k \e k )(e k \, (5) 

j>l k>l 

we order the eigenvalues as 

n > ri > ■ ■ ■ , £1 < e 2 < • ■ ■ (6) 
The minimum of Ef is then . Tj Ej , and it is reached for 

P(r)=^2r j \e j ){e j \, (7) 
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which is stationary since it commutes with H . This result 
is consistent with the extension for finite systems of the 
second law in Thomson's formulationiaLiS: the two con- 
ditions which characterize the state p(r) (commutation 
with H and ordering (JSJ) are the ones which ensure that 
this state is passive: no further work can be extracted 
from S after time r by means of cyclic processes. 

Note that, if the spectrum of H involves degeneracies 
which have no counterpart in poj the final state p(r) is 
not unique (in contrast to p cq associated with Wth)- 

Altogether, the maximum of the amount of work @ 
that can be extracted from S is 

3,k 



For W, which depends only on the initial state and 
Hamiltonian, we coin the name ergotropy ('epr/ov. work; 
rponr]': transformation, turn). By construction, we have 
Wth > W > 0. The ergotropy W vanishes if po is pas- 
sive. It equals the thermodynamical upper bound Wth 
only if there exist two numbers p and Z such that the 
eigenvalues (JSJEJ of po and H satisfy lnr, = —(3ej — \nZ, 
so as to allow p(r) = U pqW to reach a gibbsian form © 
in spite of the lack of thermalization mechanism. Notice- 
able examples include: (i) a pure initial state of S; (ii) 
two-level systems; (Hi) harmonic oscillators in case po is 
a gaussian state, since both the sequences e k and lnr^ 
are then equidistant. 

For macroscopic systems the difference Wth — W is 
typically relatively small, since the final state Q) may 
lie close to an equilibrium state if r is large; the spectra 
of In po and H are dense, and a linear relation between 
them is approximately satisfied in the relevant region. 
However, for finite systems Wth — W can be significant. 

It remains to show that the bound W can actually 
be reached by coupling S with some source of work that 
realizes a cyclic process. We thus want to find a time r 
and an interaction V{t) which vanishes at t = and at 
t = t, so that, when added to the Hamiltonian H, V(t) 
leads from the initial state po to p(r) defined in J7J). An 
evolution operator which realizes this goal is 

tf = 5>i>fol- (9) 
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A Hamiltonian H + V(t) generates in the interaction rep- 
resentation an evolution operator U\ which satisfies 

iS^ = e im ' h V(t) e- im ' h UjCt), (10) 
at 

with Ui(0) = 1. We define Ui(t) = e lHT/h U = e- lkT ' r \ 
where A is obtained by diagonalization. We choose for 
Ui{t) the simple form Ui{t) = e -^(t)/h^ where = 
ip(0) = <p{t) = and c/?(r) = r. Then the potential 
V(t) = Lp(t)e- lHt / h Ke lHt / h describes according to (0 
a source of work that extracts during the time r the work 
(jSJ) from S. This potential is far from unique. 

In case not only p{r) but also the initial state is sta- 
tionary, [po,H] = 0, we can choose the same eigenbases 
for p and H, but © implies that |r 2 ), ... in JSJ) are 
deduced from \e\), \e2), ■■■ by some permutation. Then 
the matrix | {r*j | Efc) j in JSJ is a permutation matrix. For 
instance, if the lowest two levels s± < £2 have initially 
the inverted populations r 2 < r%, respectively, so that 
\ r i) = 1^2) and |r 2 ) = |ei), we may easily implement 
the transformation U, a permutation that interchanges 
r\ and r 2 , either in a rapid or in a quasistatic regime. 

Comparison of activities. We wish to compare two 
states po and ctq of a system S as regards the maxi- 
mum work that they may provide. To make such a 
comparison meaningful, we assume the initial energies 
to be the same, E(p ) — E{ao). If S is macroscopic 
and can reach equilibrium at the end of the process, Wth 
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depends only on the entropy S of the initial state, and 
it decreases when S increases, since — dWth/dS' is the 
temperature of the final state reached. However, the sit- 
uation is different for finite systems. Consider, for in- 
stance, a three-level system with eigenenergies £1,3 = =pl, 
£2 = 0. The eigenstates of po and <jq are taken as 
|r 1)3 ) = |*i,3) = (ki> T N»A/2, |r 2 > = \s 2 ) = |s 2 ), 
so that E(pq) = E(<jq) — 0. If their eigenvalues are 
{rj} = {0.90,0.08,0.02} and { Sj } = {0.91,0.05,0.04}, 
the entropy S(p ) — 0.375 exceeds S(a ) — 0.364. Ac- 
cordingly, the thermodynamic bound W t h(po) — 0.882 
for the work is smaller than Wth(co) — 0.887. Neverthe- 
less, the ergotropy W(po) = 0.88 of p is larger than the 
ergotropy W(<tq) = 0.87 of <tq. The actually reachable 
bounds are, as expected, lower than the corresponding 
Wth's, but they are reversed in order: the entropically 
more disordered state po may provide more work. 

Thus, the entropy criterion fails for comparing the er- 
gotropies. The theory of majorization 13 , that we briefly 
recall now, provides another criterion which may be help- 
ful. In quantum statistical mechanic al 14 a density oper- 
ator p is said to majorize a if their eigenvalues r,- and Sj , 
set in the decreasing order JJJJ, satisfy 

k k 

r j > Sj, for any k > 1. (11) 

This property, denoted as p y a, is transitive (p y a 
and a y t imply p y r). It characterizes order, but 
in a stronger way than entropy since p y a implies not 
only S(p) < S(cr), but also tr/(p) < tr/(er) for any 
concave function f(x). Pure states majorize all states, 
while in a Hilbert space of dimension n all states majorize 
p = 1/n. If we have both p y a and a y p, then p and 
cr are unitarily equivalent. However, the order defined 
by majorization is incomplete since, for n > 3, pairs of 
states p and a exist of which neither majorizes the other. 

Returning to the comparison of the activities of po and 
(To with E(po) = i?((7o), we find after summation by parts 
the ergotropy difference o~W = W(po) — W(cto) as 

k 

SW = - rj )e s = 5>k+i - e k ) Yfrj - Sj ).{12) 

j>i k>i j=i 

Hence a sufficient condition for po to be more active than 
do is po y Co (S(po) < S(ao) alone is neither necessary 
nor sufficient). There exists a wide class of non- unitary 
evolutions 1 ^ which lead from a state po to Co such that 
Po >~ (Tq. If po and oo have the same energy, we have 
W(po) < W(do). For instance, the diagonal part po = 
J2k\ £ k) ( £ k\po\£k) (£k\ of po has the same energy as p 
itself, but W(po) > W(po) since p y po 13 - In order to 
find opposite behaviors in the comparison of ergotropies 
and of free energies (or entropies), we have to search for 
cases when po 7^ <jq and po )f <Jq . All pairs of inequalities 
W(po) ^ yV(co) and S(po) ^ S(cro) may then occur, as 
illustrated by examples given above and below. 



Auxiliary system. If S is supplemented with an aux- 
iliary system fi with Hamiltonian Hq and initial state 
ojq, the overall Hamiltonian H + Hq + V(t), where V(t) 
couples S + O with external sources, generates a unitary 
transformation in the product Hilbert space, which is 
more general than for work sources coupled separately 
to S and f2. The initial state po (B> wo is uncorrelated 
and the evolution conserves its (factorized) eigenvalues. 
The ergotropy satisfies W(p ® lu ) > W(po) + W(wo), 
the same inequality as for Wth- Consider, however, again 
two states po and <tq with the same energy. If they are 
macroscopic and satisfy Wth(po) > Wth (00), this order- 
ing of activities is not changed by the introduction of the 
auxiliary system, since the additivity of entropy for the 
initial states implies Wth(po ® <^o) > Wth(0o ® wo). But 
for finite systems the order of ergotropies can be reversed. 
As a first example consider for S the same three-level sys- 
tem as above, with eigenvalues of po and do now equal 
t° i r j} = {0.8,0.1,0.1} and {sj} = {0.5,0.5,0}. For 

we take a two-level system with eigenenergies and 
A > 0, initially in a pure state. Although po )/- uq, 
we have both S{p ) 0.639 < S(a ) 0.693 and 
W(po) = 0.7 > W((T ) = 0.5. However, coupling with 
f2, which does not change the entropies, reverses the in- 
equality for the ergotropies if A < 1/4, since (for A < 1) 
W(p ® uj ) - W((7 ® w ) = 0.4 A - 0.1. For A = 0, 
the auxiliary system f2 does not contribute to the energy 
balance but nevertheless raises W owing to its order. 

The opposite situation is also possible, provided S has 
at least four levels. Consider for instance a system S with 
eigenenergies £1,4 = =Fl, £2,3 = =f(1— a;), with < x < 1. 
As eigenstates of po and uo we take \ri t i) — |si,4) = 
(| £l ) ± |£ 4 ))/V2, |r 2 , 3 ) = |aa, 3 ) = (M ± \e 3 ))/V2, which 
ensures E(po) = E(<jq) — 0, and as eigenvalues: 

{ri} = ^ {^ + 3)^,^,0}, (13) 

= ITT \2 i 2w > 2w ' w ^ " w ) 2 < (! - w ) 3 } » ( 14 ) 
(1 + wy 

which are ordered according to © provided 1 > w > 1/2. 
The fact that po )/- ctq and (To )/- po allows to violate the 
thermodynamical ordering, since we have simultaneously 
S(po) < S(a ) and W(p ) < W(ct ) for sufficiently small 
x, as seen from (l + w) 2 [W(a ) -W(p )] = (l-w)(2w- 
1) — xw(l + 2w — w 2 ). Take as a two- level system; its 
only relevant feature will be the eigenvalues {w, 1 — w} 
of its initial state ujq. Provided 1/4 < w 2 < 1/2, the 
equations (f 1 1|) for k = 1,2,..., 8 are satisfied and hence 
Po®ojo y oq®ujq. This implies W(po<8> o>o) > W(<To®wo): 
n restores for the ergotropies the order inferred from the 
thermodynamic relation S(po) < S(ao)- 

All these contradictions between the predictions of 
thermodynamics and the behavior of finite systems do 
not occur for the subset of states that may ordered in 
the sense of majorization: then p y cr implies both 
S(j>o) < S{a a ) and W(p ) > W(cr ) and W(p ® u ) > 
W((To <8> wq) for an arbitrary loq. 
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Correlations. In the above examples the initial state 
of S+O was uncorrelated, but not the final state, because 
evolution permutes its eigenvectors although its eigenval- 
ues remain factorized. Thus, contrary to thermodynamic 
intuition, the maximum work can be achieved owing to 
creation of correlations. 

Conversely, if the initial state Qo of S+f2 is correlated, 
thermodynamics predicts that the available work is in- 
creased, due to the subadditivity of entropy, S(Qo) < 
S(po) + S(cto), where p = tr f2 Qo and cr = tr s Qo are 
the marginal states of S and 17, respectively. However, 
for finite quantum systems, we have to compare the er- 
gotropies W(Qo) and W(po (8 0o). Take for S a system 
with three energy levels £j, i = 1,2,3 and for O a two- 
level system k = 1,2 with energy levels 0,£ such that 
< e < £2 — £i, e < £3 — £2, and for S + SI a station- 
ary initial state Qo- Denoting the common eigenstates 
of H + Hq and Qo as \i,k), we assume that the eigen- 
values qik of Qo are ordered as qn > q\ 2 > q 2 \ > q 22 > 
931 > <?32- Then Qo is both correlated and passive. Sup- 
pressing correlations leads to a factorized state po <8> 0J0 
with the same energy, and eigenvalues = Y^t=i Qik 
for po and pk — X)i=i Qik f° r w o> The ordering of the 
set riPk may now differ from that of qik', for instance, 
if <7n is close to one and all other qik's are small with 
the same order of magnitude, we have r2P2 < T^pi, and 
W(pq<Z>ujo) = (s3-e 2 -e)(r 3 p 1 -r 2 p2) > 0. Suppressing 
the correlations has thus let the system S + 11 become 
active. Altogether, the order carried by correlations, al- 
though it manifests itself directly in the entropy, has no 
relation with increase or decrease of ergotropy. 

Fluctuations of work. The work —W extracted from S 
by means of external sources acting in a cycle was defined 
by as the expectation value of the difference between 
the initial and final energies of S. This expectation value 
can also be observed as an averageii by coupling the 
macroscopic sources of work to N ^> 1 identical systems 
S operating in parallel, in such a way that their total 
Hamiltonian is the sum of the Hamiltonians H + V (t) of 
each one. However, for a single system S, the work is 
expected to fluctuate. In order to evaluate such fluctua- 
tions, we need to define work as a random variable for the 
dynamics described by the optimal unitary transforma- 
tion ©. To do this, we need to determine, for each single 
realization in the statistical ensemble described by po, the 
initial and final energy of S. This can be done unambigu- 
ously, without perturbing the overall ensemble po, only if 
both po and p(r) commute with H. Then the measure- 
ment of energy at the time t = allows to filter systems 
which have the same initial energy si' after the mea- 
surement. The evolution operator J^J then determines 
for each initial £?' the final energy £, , hence the work 
—e^p is given by S to the sources. We can thus regard 



£? — e^p as a classical random variable describing work, 
and governed by the diagonal elements rj of the density 
matrix p . In the same way as W = J2jk r jPjk(£k — 
where Pjk = \ (sj\U\ek)\ 2 is a permutation matrix, we get 
the ergotropy dispersion: 

AW 2 = J2 rjP jk (e k - e 3 ) 2 - W 2 . (15) 
jk 

It vanishes for passive states, for which Pjk — Sjk and 
W = max W = 0: one cannot hope to put no work 
in average, and get some work out by making use of 
fluctuations. For active states, AW may become size- 
able. In the example of a two level-system with the ini- 
tial state: p = ri|e 2 )(£2| + ?*2|£i)(£i|, the ergotropy is 
W = (ri — r 2 )(£2 — £1) and its fluctuation from l(l5|) is 
AW= 2^/ryr2~(e 2 - £1). The relative fluctuation AW/W 
vanishes for r\ — 1, r 2 =0, when W is largest; it be- 
comes large when W is small, i.e., r\ — » r 2 — > \. Many 
systems S should then operate in parallel to produce on 
average the amount W in a predictable fashion. 

Conclusion. Maximal work extraction is one of the 
basic problems of thermodynamics and has applications 
in various processes of energy conversioni2i&. In macro- 
scopic physics, the answer (0J is governed by the non- 
decrease of entropy. We have shown that finite devices 
are less efficient in this respect: any evolution of a ther- 
mally (but not mechanically) isolated quantum system 
must leave unchanged not only this entropy, but all the 
eigenvalues of the density operator, which prevents in 
most situations the thermodynamic bound from being 
attainable. We have given a general explicit expression 
JSJ) for the "ergotropy" , the quantity which, for finite sys- 
tems, replaces the free energy: it is the upper bound of 
the work that a finite system S in a non-equilibrium ini- 
tial state po may yield if it is coupled to external sources 
of work undergoing a cyclic transformation. Many inter- 
action Hamiltonians V(t) allow to reach this bound. 

The proper measure of order for comparing the abili- 
ties of work production of finite systems is thus ergotropy 
and not free energy. Several consequences of this result 
contradict thermodynamic intuition. Consider, for in- 
stance, a state o~o of S having the same energy as po and 
lower entropy. Thermodynamics suggests that more work 
might be extracted from do than from po, and moreover 
that the presence of an auxiliary system 17, in a state 
wo initially uncorrelated with S, preserves this property. 
Such statements can be violated in finite quantum sys- 
tems. However, even for finite systems, there is a domain 
of states defined by the majorization relation, where pre- 
dictions of thermodynamics are qualitatively correct. 
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